Unusual conductance of polyyne-based molecular wires 
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We report a full self-consistent ab initio calculation of the current-voltage curve and the con- 
ductance of thiolate capped polyynes in contact with gold electrodes. We find the conductance of 
polyynes an order of magnitude larger compared with other conjugated oligomers. The reason lies 
in the position of the Fermi level deep in the HOMO related resonance. With the conductance 
weakly dependent on the applied bias and almost independent of the length of the molecular chain, 
polyynes appear as nearly perfect molecular wires. 
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The study of transport properties of single molecules 
have attracted a significant attention because of their 
potential use in molecular electronic devices. One 
of the major classes of molecules considered in con- 
ductivity studies, primarily for their molecular wire 
behaviorfi*2iiii£ is conjugated oligomers. They have 
shown a number of useful nonlinear properties such 
as conductance switching and negative differential 
resistances^ However, inspite of a number of interest- 
ing experimental a molecule with good molecular wire 
properties has not yet been spotted. 

A useful molecular wire should provide a high and sta- 
ble conductance over a wide bias region and for various 
lengths of the molecules. A linear chain of carbon atoms 
with double bonds between neighboring atoms, usually 
referred to as cumulene, was proposed as an ideal molec- 
ular wirei and the calculations of the conductance of cu- 
mulene connected to gold electrodes were reported. 9 ' 10 
Lang and Avouris 9 - showed that the conductance of cu- 
mulene did not stay constant in the ballistic regime, but 
rather oscillated between the constant values character- 
istics of the odd and even number of atoms in the chain. 

In this paper we show an entirely different behavior 
of polyynes, another form of the carbon atom chain. 
Polyynes are simple and yet most intriguing of conju- 
gated organic oligomers. Only recently, have they been 
assembled up to decayne^i Formed as a linear chain of 
carbon atom pairs (CC)„, with alternating single and 
triple bonds, they are a unique, truly one-dimensional, 
molecular system. Two 7r-electron systems of the sp- 
hybridized structure provide polyyne with approximately 
cylindrical electronic derealization along the conjugated 
backbone. The electronic transport is therefore indepen- 
dent of the rotation around the single bond, which is a 
limitation often present in other organic oligomers. 7 

We have obtained the electronic structure and trans- 
port properties of a series of polyynes up to octayne, 
connected to gold electrodes. The stability of polyyne 
with respect to single- and triple-bond alternations was 
achieved by fixing the molecule at the ends with thiol 
bonds. In addition, the thiol capped polyynes make a 
strong chemisorption bond onto the metallic electrodes. 
We found that they had more than an order of mag- 



nitude higher conductance when compared with other 
conjugated oligomers. In contrast to the cumulenes they 
were not prone to oscillations in conductance with the 
length of the molecule. We also found that their conduc- 
tance was very weakly dependent on the molecular length 
and almost constant in the wide range of bias voltages. 

A previous study, with a different theoretical ap- 
proach, included polyynes up to triyne with Pd contacts 
of a different contact geometry which resulted in lower 
conductance^ 

In order to perform the first-principle quantum model- 
ing of the electronic structure under nonequilibrium con- 
ditions and to calculate the current-voltage characteris- 
tics and the differential conductance of the system, a full 
self-consistent ab initio method, which includes portions 
of the electrodes, had to be usedii^ The calculations were 
carried out using a nonequilibrium Green functions tech- 
nique based on density functional theory, as implemented 
in the TranSIESTA simulation packaged The current 
through the contact was calculated using the Landauer- 
Buttiker formula, 15 I(V b ) = G T{E, V b ) dE, where 

Go = 2e 2 /h is the quantum unit of conductance and 
T(E, Vb) is the transmission probability for electrons in- 
cident at an energy E through the device under the poten- 
tial bias V>. The difference between the electrochemical 
potentials Hl/r, of the left and right electrodes, respec- 
tively, is ~ — eVb- The computational procedure 
used was described extensively elsewhere^. 

The molecular electronic system considered consisted 
of a monolayer of molecules coupled to two semi-infinite 
electrodes, as depicted in Fig^ for hexayne. We opti- 
mized the geometry of free thiol capped polyyne in a 
separate DFT calculation 1 ^. The molecule was then po- 
sitioned perpendicularly to the z-direction in the hollow 
sites of both Au(lll) surfaces of the electrodes symmet- 
rically at a favorable Au-S bonding distance^. We let 
the molecular coordinate relax, keeping the gold atoms 
at their bulk positions. A small change in the molecular 
geometry occurred with respect to the free molecule. We 
discuss this at the end of the paper. 

The main characteristics of the obtained transmission 
spectra are the almost linear increase of the current with 
bias voltage and a high value of conductance over a wide 




tmm 




FIG. 1: Hexayne connected to two Au(lll) surfaces via thi- 
olate bonds, shown with the MPSH HOMO orbital (upper 
panel, notice the delocalized shape at each side of single bond) 
and the MPSH LUMO orbital (lower panel) at V b = 0.6V bias 
voltage. 
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FIG. 2: Current and differential conductance of polyyne sys- 
tems as a function of bias voltage. A self-consistent calcula- 
tion has been employed for each bias voltage. Notice the slow 
and smooth almost linear decrease of the conductance. 



bias region from -2V to 2V, as seen in Fig|2j The ob- 
tained spectra of all members^ of the polyyne series show 
a mutual pronounced similarity without oscillations with 
the length of the molecule. The conductance at zero bias 
voltage is 1.65 Go for diyne, 1.56 Go for tetrayne, 1.49 
Go for hexayne, and 1.44 Go for octayne. The overall 
high value of conductance is expected, since the main 
transmission channels involve double degenerate molec- 
ular 7r orbitals. Unlike in the true ballistic transport 
in a quantum structure where a constant conductance 
might be expected irrespective of the molecular length^, 
a slow decrease of conductance with the length of the 
molecule is demonstrated. The decrease should be at- 
tributed to a weak reduction of hybridization at the 
molecule-electrode contact with the number of atoms in 
the molecular chainSS. 

When compared with other conjugated oligomers, 
polyynes show much higher conductance. In Fig|3| 
we compare the transmission amplitude at zero bias 
of three different thiolate capped molecules: hexayne, 
diphenyl diacetylene (DPA2), and phenylene vinylenc 
oligomer with three benzene rings (OPV3). They were 
all chemisorbed onto Au(lll) electrodes in hollow posi- 
tions at both ends in order to ensure the same bonding 
geometry thus avoiding the possible bonding site effects 
on the considered transmission. The resulting interelec- 
trode separations for all three systems were rather simi- 
lar, but the corresponding zero bias conductances differed 
by more than an order of magnitude, as shown in Table 
[U The large difference in conductance is a consequence of 
the entirely different electronic structure and density of 
states at the Fermi level (Ef). The hybridization of the 
molecular level with the gold electrode states broadens 
the level into the resonance. The width and the position 
of the resonance with respect to the Fermi level of the 
system depend on the internal structure of the molecule 



and its bonding to the electrode. In the hexayne case, the 
double degenerate molecular tt orbitals are involved. The 
resulting HOMO related resonance 21 is wide and shifted 
up in energy close to the Fermi level, as seen in Fig|3] In 
fact, the HOMO and LUMO related resonances strongly 
overlap giving a wide transmission band with large den- 
sity of states, which results in high transmission. The 
transmission of DPA2 and OPV3 systems is consider- 
ably smaller when compared with the hexayne, as shown 
in FigEl The reason for that is the lifted HOMO level 
degeneracy in DPA2 and OPV3 and the positioning of 
the HOMO level further below Ef, which results in the 
lower density of states at the Fermi level. 

The high density of states at the Fermi level is the 
property of other members of the polyyne family, as well. 
With the increased length of the polyyne chain, levels are 
becoming more densely distributed and of shorter width. 
The remarkable property is, however, that the position 
of the HOMO related resonance moves slightly closer to 
Ef, thus compensating for the decrease of the density of 
states at Ef due to the level sharpness, as seen in FigQ] 
That results in a very weak dependence of the zero bias 
conductance on the molecular length. 

The conductance of polyynes decreases smoothly with 
the increase of the bias voltage, as shown in Fig|2 In or- 



TABLE I: Interelectrode separation d e i ec and zero-bias con- 
ductance G([i,0) of polyyne systems compared with molecules 
of similar length. d e i ec is the distance between the surface gold 
planes of the electrodes. 
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FIG. 3: Transmission amplitude of hexayne compared with 
DPA2 and OPV3 at zero bias voltage. MPSH eigenvalues 
are marked with circles for hexayne, diamonds for DPA2, and 
triangles for 0PV3. 



der to elucidate this behavior, we show the transmission 
amplitudes of hexayne for a set of bias values in Fig|S] A 
remarkable similarity, in both shape and value, of trans- 
mission followed by a slight shift of the position of LUMO 
and higher resonances with the bias is seen. The sepa- 
ration between HOMO and LUMO levels decreases only 
by 0.11 eV at the bias of 1.8 V with respect to the zero 
bias value. The consequence of the similarity of trans- 
mission spectra is the relatively weak dependence of the 
conductance of hexayne on the bias, as seen in Fig|2j 

The other polyynes showed the dependence of trans- 
mission amplitudes when subjected to different bias volt- 
ages similar to those for hexayne. With the change of the 
bias from V to 2.0 V, the HOMO-LUMO energy sepa- 
ration for all polyynes that we consider decreases by no 
more than 0.15 eV. Such a small change in the molec- 
ular level positions relative the average electrochemical 
potential of the electrodes ensures a slow decrease of the 
conductance over the entire bias region, as clearly seen 
in FigUl 

What makes polyynes so much different from other 
conjugated polymeres is the high density of states at 
the Fermi level. Even in a free molecule of polyyne, 
where the alternate single and triple bonds open up 
a HOMO-LUMO gap, the strong polarizability and 
hyperpolarizabilityliSS evolve, owing to the large num- 
ber of electrons. When thiol capped polyynes are 
chemisorbed onto gold electrodes, the strong hybridiza- 
tion of molecular states with the metallic electrode states 
results in a new electronic structure of the combined sys- 
tem. The change of the level positions is so strong that 
the Fermi level of the system enters deep into the HOMO 
related resonance resulting in the pronounced metallic 
character of the system. 

As already said, we used LDA for the exchange- 
correlation functional in our DFT calculations. There 
has been a discussion in the literature on the validity of 
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FIG. 4: Transmission amplitude of polyynes at zero bias 
voltage. Note the rather stable position of HOMO related 
peaks for all polyynes. 

DFT with the local density (LDA) and generalized gradi- 
ent (GGA) approximations for 7r-conjugated systems^. 
Here, however, the 7r-conjugated system is connected to 
the electrodes and shows more metalliclike character. 
LDA- and GGA-based DFT calculations are therefore 
expected to give a respectable accuracy usually obtained 
in DFT calculations. The crucial parameters determin- 
ing the HOMO-LUMO gap are the carbon-carbon bond 
lengths and their alternations. We performed the calcu- 
lations for isolated molecules of thiolate capped polyynes 
using the Gaussian03 program^ with several XC func- 
tionals and found that the GGA results did not substan- 
tially improve the values of the gap, nor the single and 
triple bond lengths. Indeed, the difference in the HOMO- 
LUMO gaps was less than 0.037 eV between the PW91 
(GGA) and SVWN (LDA) based calculations, as seen in 
Table ITT1 When extrapolated to infinite chain length the 
gap stayed open, in agreement with the calculations with 
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FIG. 5: Transmission amplitude of dithiolhexayne connected 
to gold electrodes as a function of bias voltage. Energies mea- 
sured from the average electrochemical potential of the elec- 
trodes. 
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TABLE II: MPSH HOMO-LUMO energy gap AE, in eV, 
at zero bias of thiolate capped polyynes with Au electrodes 
compared with the HOMO-LUMO gaps of bare thiol capped 
polyynes (DFT calculations with SVWN, PW91, and PBE 
XC functionals) . 

AuS{CC) n SAu HS(CC) n SH 

n AE Ehomo Elumo AEsvwn AEpw9i AEpbe 



2 


3.254 


-0.897 


2.357 


3.529 


3.559 


3.549 


4 


2.063 


-0.643 


1.420 


2.333 


2.370 


2.364 


6 


1.520 


-0.482 


1.038 


1.762 


1.797 


1.788 


8 


1.211 


-0.385 


0.825 


1.428 


1.463 


1.458 



more exact exchange potentials of Weimer et ali^ For the 
molecule bonded to the electrodes, the HOMO-LUMO 
MPSH gap (first column of Table [Hj was reduced with 
respect to the isolated molecule as expected. The reduc- 
tion was slightly more pronounced for shorter molecules, 
equal to 0.031 eV for diyne, and decreased to 0.025 eV 
for octayne. Evidently, the HOMO-LUMO MPSH gap 
closely followed the gap of the isolated molecule. The 
possible underestimation of the gap in the LDA-based 
calculations compared with the nonlocal XC potential 
calculations would have no major effect. The metalli- 
clike character of the conductance appears to be the con- 
sequence of the hybridization of the 7r-electron systems 
of polyynes with the gold electrode states, irrespective of 
the details of the XC potential. 

The hybridization with the gold states did not affect 
the bond alternation in the electrode-connected hexayne, 



as seen in Table ITTT1 The slight reduction of single bonds 
and the expansion of triple bonds of the connected hex- 
ayne are noticed. The change of the bond length is the 
interplay of the hybridization of the molecular states with 
the gold electrode states and the relaxation of position of 
molecular atoms. The amount of change is smaller when 
the molecule is longer owing to the redistribution of the 
relaxation over the entire molecule. 



TABLE III: Average single and triple related bond distances, 
measured in A, of thiolate capped hexayne in contact with Au 
electrodes compared with the bare thiol capped hexayne bond 
distances. 





AuS(CC) 6 SAu 


HS(CC) 6 SH 








SVWN PW91 


PBE 


^single (-^-) 


1.316 


1.321 1.330 


1.332 


dtriple (A) 


1.286 


1.249 1.252 


1.253 



In this paper, we have shown that polyynes are 
molecules with high conductance over the large bias re- 
gion almost irrespective of the length of the molecule con- 
sidered. Although polyynes form 7r-conjugated systems, 
they show an almost metalliclike character of the current 
transport, when connected to the electrodes. This makes 
them a possible candidate for a good molecular wire to 
be considered in molecular nanoclcctronics. 
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